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We present new measurements of the six largest branching fractions of the Ki^ using 
data collected in 1997 by the KTeV experiment (E832) at Fermilab. The results are 
B{Kl -^ 7r±e=Fi/) = 0.4067 ± 0.0011, B{Kl -^ n^/i^u) = 0.2701 ± 0.0009, B{Kl -^ -k+tv-tv°) = 
0.1252±0.0007, B{Kl -^ tt^W") = 0.1945±0.0018, B{Kl -> tt+tt") = (1.975±0.012) x 10"^ and 
B{Kl -^ tt^tt") — (0.865 ±0.010) x 10~^, where statistical and systematic errors have been summed 

in quadrature. We also determine the CP violation parameter \r}^ | to be (2.228 ± 0.010) x 10""^. 

Several of these results are not in good agreement with averages of previous measurements. 
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I. INTRODUCTION 



Most recent experimental work on the K^ has focused 
on rare and CP violating decays. The branching frac- 
tions of the main K^ decay modes, however, have not 
been measured together in a modern, high-statistics ex- 
periment. These branching fractions are fundamental ex- 
perimental parameters used to determine the GKM ele- 
ment \Vus\^ the CP violation parameter \'q^ |, and the 

normalization for many other rare decay measurements. 

In this paper, we present new results for the six largest 



measuring the following ratios of decay rates: 



Kl branching fractions: Ki^ 
' Kl -^ ttOttO 
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We determine these branching fractions by 
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Each ratio is measured in a statistically independent 
data sample collected by the KTeV (E832) experiment 
at Fermilab. Note that throughout this paper, inner 
bremsstrahlung contributions are included for all decay 
modes with charged particles. 

Since the six decay modes listed above account for 
more than 99.9% of the total decay rate, the five partial 
width ratios may be converted into branching fraction 
measurements. For example, the Kj^ -^ tt e^v branch- 



ing fraction, BkcS, may be written as 
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(6) 

where B^are = 0.07% is the sum of branching fractions 
of other rare Kl decay modes [i^- In terms of our 
measured partial width ratios, T^^/Txes = Too/rooo x 

rooo/rj<'e3- 

The paper is organized as follows. First, we give a brief 
description of the KTeV detector and the data sets used 
in this analysis. Next, we present an overview of the 
analysis techniques followed by a more detailed discus- 
sion of selection criteria for the individual decay modes. 
Section IVIII contains a description of the Monte Carlo 
simulation used to determine the detector acceptance, 
and a discussion of the resulting systematic uncertain- 
ties. In Section rVlllI we present results for the partial 
width ratios and branching fractions, along with several 
crosschecks of our analysis. Our branching fraction mea- 
surements are then used to extract \ri^ |. Finally, our 

results are compared with previous measurements. 



II. KTEV DETECTOR 

The KTeV detector (see Fig. ^ and associated event 
reconstruction techniques have been described in detail 
elsewhere [l| . Here we give a brief summary of the essen- 
tial detector components. An 800 GeV/c proton beam 
striking a BeO target is used to produce two almost par- 
allel neutral beams, shaped by a series of collimators. 
Except for a special "low-intensity" run described later, 
a fully-active regenerator is placed in one of the beams 
to provide a source of Ks for the measurement of e'/e; 
decays in the "regenerator" beam are not used in this 
analysis. The other beam, referred to as the vacuum 
beam, provides K^ decays used for these measurements. 
A large vacuum decay region extends to 159 m from the 
primary target. 

Following a thin vacuum window at the end of the 
vacuum region is a drift chamber spectrometer used to 
measure the momentum of charged particles; this spec- 
trometer consists of four chambers, two upstream (DCl- 
DC2) and two downstream (DC3-DC4) of an analysis 
magnet that imparts a 0.41 GeV/c momentum kick in 
the horizontal plane. Each chamber measures positions 
in both the x and y views (transverse to the beam di- 
rection). Farther downstream lies a trigger (scintillator) 
hodoscope to identify charged particles, and a 3100 crys- 
tal, pure cesium iodide (Csl) electromagnetic calorime- 
ter. Downstream of the Csl calorimeter there is a muon 
system consisting of scintillator hodoscopes behind 4 m 
and 5 m of steel. Veto detectors surround the vacuum 
decay region, each drift chamber, and the Csl calorimeter 
(Vacuum- Veto, Spec- Veto, and Csl- Veto). 

KTeV uses a three-level trigger system to reduce the 
total rate of recorded events. The Level 1 and Level 2 



triggers are implemented in hardware and the Level 3 
trigger is a software filter that uses the full event recon- 
struction. With the exception of Kl -^ tt'^tt'^, the data 
samples used in this analysis do not require the Level 3 
trigger. 

The analysis of the Kl branching fractions benefits 
from the KTeV detector design that was optimized to 
measure the direct CP violation parameter using a Monte 
Carlo simulation (MC) to determine the detector accep- 
tance. To reduce uncertainties in the simulation, it is 
important that apertures and detector geometry are well 
measured, and that there is very little material before 
the calorimeter to affect decay products. 

The Z-locations of detector elements are known with 
100 /im precision. The transverse sizes of the drift cham- 
bers are known to 20 /im; the transverse dimensions and 
relative locations of the other detector elements are de- 
termined to better than 200 /im using charged particle 
tracks. 

The detector has very little material upstream of the 
calorimeter, reducing losses from multiple scattering, 
hadronic interactions, and 7 — > e+e^ conversions. The 
material from the vacuum window to the last drift cham- 
ber is only 0.012 radiation lengths (Xq), or 0.007 pion 
interaction length (Aq); the material downstream of the 
last drift chamber to the front face of the Csl is 0.031 ATq, 
or 0.014 Aq. 

The analysis presented in this paper also benefits from 
extensive detector calibration performed for the e'/e anal- 
ysis [l| . The momentum kick of the analysis magnet is 
determined to 0.01% precision using Kl — * tt+tt^ events 
and the PDG value of the kaon mass j2|. The Csl 
calorimeter energy scale is determined to better than 
0.1% based on calibration using 500 million momentum 
analyzed electrons from Kl —^ ir^e^i/ events. The mo- 
mentum resolution of the spectrometer and the electro- 
magnetic energy resolution of the Csl calorimeter are 
both better than 1%. 



III. DATA COLLECTION 

In this paper, we report results based on data taken 
during two periods in 1997, which will be referred to 
as "high intensity" and "low intensity." All five partial 
width ratios are measured in both the high and low in- 
tensity data samples. The quoted result for each ratio 
is based on the sample yielding the smaller total uncer- 
tainty; the other sample is used as a crosscheck. With 
the exception of TKus/^Ke Si t he high intensity sample 
gives a smaller uncertainty |4lj . 

The high intensity period was used primarily to col- 
lect AT ^ TTTT decays for the i?e(e'/e) measurement ij. 
In addition to the nominal 2-pion triggers, several addi- 
tional triggers with relaxed requirements were included 
for systematic studies; these additional triggers provide 
the data samples for the branching fraction analysis^3- 

The low intensity data were collected during a 2-day 
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FIG. 1: Plan view of the KTeV (E832) detector. The evacuated decay volume ends with a thin vacuum window at Z = 159 m. 
The label "Csl" indicates the electromagnetic calorimeter. 



special run. The beam intensity was lowered by a fac- 
tor of 10, and the drift chamber operating voltage was 
raised to increase efficiency. In addition, the regenerator 
was removed, resulting in two vacuum beams and elim- 
inating extra detector hits from the interaction of beam 
particles in the regenerator. The data collected during 
the low intensity period have significantly lower detector 
activity. For example, the average number of spurious 
drift chamber hits is only 2.3 in the low intensity period 
compared to 43 for the high intensity period. 



IV. MEASUREMENT STRATEGY 



requirements may be used to distinguish different kaon 
decay modes from each other, and to reduce background 
to a negligible level for all decay modes. For some de- 
cay modes, the excellent spectrometer and calorimeter 
resolution allows us to achieve this background rejection 
without using all of the available detector information. 



For the K^ -^ n^fi^i^ and K^ 



IT decay modes. 



we exploit this flexibility to reduce systematic uncertain- 
ties in the acceptance. 

For the K^ — > tt jjl^v decay mode, we do not make use 
of the muon system to identify the muon; this avoids sys- 
tematic errors in modeling muon propagation in the steel 
in front of the muon hodoscope as well as in modeling 



The analysis is optimized to reduce systematic uncer- 
tainties resulting from the Monte Carlo simulation used 
to correct for acceptance differences between pairs of de- 
cay modes. With the exception of Fooo/FKes, we consider 
ratios of decay modes with similar final state particles. 
To be insensitive to the absolute trigger efficiency, events 
for the numerator and denominator of each partial width 
ratio are collected with a single trigger (the only excep- 
tion is the measurement of Foo/Fqoo 43]). 

The trigger requirements for each partial width ratio 
are summarized in Tabled The main trigger requirement 
for each pair of decay modes is either two charged tracks 
or a large energy deposit in the Csl calorimeter. For 
rooo/rife3) the trigger requirement is only 25 GeV of 
energy in the calorimeter for both decay modes; there is 
no charged-track requirement for K^ — *• n e^v. 

As will be described below, very simple event selection 



TABLE I; Trigger requirements used to measure each par- 
tial width ratio. Note that two different triggers are used to 
measure roo/Tooo. "Two charged tracks" refers to hits in the 
drift chambers and trigger hodoscope, "total Csl energy" re- 
quires more than 25 GeV energy-sum in the 3100 channels, 
"Csl clusters" refers to the number of clusters above 1 GeV, 
and "vetos" are used to reject events. 
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gaps between the muon counters. For the Kl — > 7r+7r~7r° 
decay mode, tt'^ — > 77 is not reconstructed in order to 
be insensitive to pion showers that bias the photon en- 
ergy measurement in the Csl calorimeter. Also, by ig- 
noring the tt", the reconstruction of Kl —> TT~^TT~Tr^ and 
Kl -^ IT e^v decays are very similar, reducing the un- 
certainty on the acceptance ratio of these two modes. 
As a crosscheck (Section IVIIIBp . the Kl — > ^T^^JL^v and 
Kl — > TT+TT^Tr" modes are also analyzed using the muon 
system and fully reconstructing the tt" . 

All reconstructed decay modes are required to have 
kaon energy, Ek, between 40 and 120 GeV, and decay po- 
sition, Zx, between 123 and 158 m from the target. For 
the reconstruction of semileptonic and Kl -^ vr+Tr^vr" 
decays, there is a missing particle {v or tt*^); this leads 
to multiple kaon energy solutions. All energy solutions 
are required to be in the accepted range. These Ek and 
Zk ranges are more restrictive than in the e'/e analysis 
in order to have a negligible contribution from Ks -^ tttt 
decays. 

In the following two sections, we discuss the techniques 
used to reconstruct "charged" decay modes with two 
oppositely-charged particles, and "neutral" decay modes 
with only photons in the final state. More details of the 
KTeV event reconstruction are given in [l|. 



V. CHARGED DECAY MODE ANALYSIS 



Charged Decay Mode Reconstruction and 
Event Selection 



The reconstruction of Kl 

^0 



TT e 



T,. 



r±/,.=F, 



i/, Kl -^ -K^ii^v, 
Kl -^ vr+vr^Tr", and Kl — > tt+tt^ begins with the iden- 
tification of two oppositely charged tracks coming from 
a single vertex. To pass the event selection, one of the 
two tracks must be within 7 cm of a Csl cluster; the sec- 
ond track is not required to have a cluster match. As will 
be discussed later, this relaxed track-cluster matching re- 
quirement reduces the inefficiency arising from hadronic 
interactions upstream of the calorimeter. 

The transverse (A, Y) decay vertex is required to 
be within an 11x11 cm^ square centered on the beam 
(RING < 121 cm^ lill). The beam profile is about 10x10 
cm^ at the Csl calorimeter. This cut removes most events 
in which the kaon has scattered in a collimator. 

The fiducial region for the charged decay modes is de- 
fined by requiring that projections of tracks fall safely 
within the boundaries of the drift chambers, trigger ho- 
doscope, and the Csl calorimeter. The tracks are also 
subject to a "cell separation" cut 1] , which requires that 
the tracks never share the same drift chamber cell. This 
requirement introduces an effective inner aperture, re- 
jecting pairs of tracks with a very small opening an- 
gle. At the Csl calorimeter, the tracks are required to 
have a large 40 cm separation to minimize the overlap of 
hadronic and electromagnetic showers. 

The different charged decay modes are distinguished 



from each other on the basis of particle identification and 
kinematics. The calorimeter energy measurement {E), 
combined with the spectrometer momentum (p), is used 
to distinguish electrons and pions. Figures [2Ia) and (b) 
show data and Monte Carlo (MC) E/p distributions for 
electrons and pions, respectively. Electron candidates are 
required to have E/p greater than 0.92; this cut retains 
99.8% of the electrons and rejects 99.5% of the pions. 
Pions are required to have E/p less than 0.92. 

In the Kl — *■ tt^ n^v analysis, the E/p requirement is 
used to reject ATeS decays. We also require that at least 
one track point to a Csl cluster with energy less than 
2 GeV (Fig. Hfc)). The 2 GeV cluster requirement re- 
tains 99.7% of the A'^a signal. This requirement does 
not distinguish the pion from the muon in 1/3 of the 
events, because 1/3 of the pions do not shower in the Csl 
calorimeter. Since the pion and muon are not identified, 
there are four Ek solutions that are all required to be 
within the 40-120 GeV range. Recall that muons are not 
identified with the muon system in order to reduce accep- 
tance uncertainties; for the other decay mode analyses, 
however, the muon system is used in veto to suppress 
background from AT^a and pion decays. 
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FIG. 2: (a) E/p (calorimeter energy over spectrometer mo- 
mentum) distribution for electrons; (b) E/p distribution for 
pions; (c) minimum cluster energy deposit among two tracks 
for A^3 decays (events with energy deposits greater than 
5 GeV are shown in the last bin). The arrows indicate se- 
lection requirements. 

In addition to particle identification for the two tracks, 
kinematic requirements are also used to distinguish the 
charged decay modes. We use three variables (m^Trj 

pI, and fc_| o)j each computed under the assumption 

that both tracks are charged pions. The only fully re- 



constructed decay, K^ ^ Tr~^Tr~ , is isolated using the 
two-track invariant mass (wi^jr) and using the two-track 
transverse momentum-squared (p^) measured with re- 
spect to the hne connecting the primary target and decay 
vertex. To separate K^ — > tt+tt^tt'^ from semileptonic 
decays, we use an additional kinematic variable. 
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where rriK and 771^^0 are the kaon and n^ masses, respec- 
tively. For Kl -^ TT+vr^vr'' decays, k^ corresponds to 

the square of the longitudinal momentum of the tt" in 
the reference frame in which the sum of the charged pion 
momenta is orthogonal to the kaon momentum. 

To illustrate the use of these 3 variables, Fig. O shows 
data and MC distributions of tti^tt, Pt, and fcn — for 
all two-track events before particle identification and 
kinematic requirements are applied. The different MC 
samples are normalized to each other using the branch- 
ing fractions reported in this paper. Figures Ofa) 
and (b) show peaked distributions in ttItttt and p^ for 
Kl -^ TT+TT^ decays. K^ -^ tt+tt^ candidates are se- 
lected with 0.488 < m^^ < 0.508 GeV/c^ and p'^ < 
2.5 X 10~^ GeV^/c^; the combined efficiency of these 
two requirements is 99.2%. For the other charged de- 
cay modes, the same 771.^,^ requirement is used to veto 
background from K^ — > tt+tt" decays. 

Figure |3fc) shows the fc_| distribution. There 



correspond- 
^0 



are two well-separated enhancements, 
ing to semileptonic and Kl -^ 7r"'"7r~7r'^ decays. 
We require K^ -^ tt~^tt~tt'^ candidates to have 
fc+_o > -0.005 GeVVc^ which retains 99.9% of 
the Kl -^ vr+Tr^vr" decays. For semileptonic decays, we 
require fc+_o < -0.006 GeVVc^; this keeps 99.5% of 
Ke3 decays and 97.9% of K^^ decays. 

In addition to misidentifying kaon decay modes, back- 
ground can also arise from hyperon decays. Back- 
ground from A — > p7r~ decays is suppressed to a neg- 
ligible level for all charged decay modes by removing 
events with invariant p'K~ mass consistent with toa 
(1.112 — 1.120 GeV/c^); to determine nipTr, the higher 
momentum track is assumed to be the proton. 

Figures BlBl show the data and MC distribution of niTrTr , 
Pf, and k^ for each charged decay mode, after apply- 
ing all event selection requirements except for the re- 
quirement on the plotted variable. The data and MC 
distributions agree well; the background for each signal 
mode is based on simulating the other three charged de- 
cay modes. 

In addition to the general reconstruction and event se- 
lection described above, several notable features specific 
to certain decay modes are described below. 

The Kl -^ n^e^v decay mode is used in four of the 
five partial width ratios. Some Kl — * n e^iy selection 
requirements are adjusted depending on the ratio. For 
^Kfj.3/^Ke3, the pion is required to satisfy a stricter E/p 
requirement {E/p < 0.85 instead of 0.92) to reduce back- 
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FIG. 3: For all two-track events, without particle identifica- 
tion or kinematic requirements, distributions are shown for (a) 
m^„, (b) Pf, and (c) k^ o- Data are shown as dots. MC simu- 
lations for Kl — >■ vr e^v, Kl —^ tv IJ^v, Kl — > 7r^7r~7r , and 
Kl — + tt^tt" are indicated. The sum of all MC contributions 
is shown as a solid line. 



ground from Kf.3 in the K^^ sample. This stricter E/p 
requirement retains 99.1% of pions, and rejects 99.93% of 
the electrons. Since there is one pion in the final state for 
both K^3 and Ke3, we use the same strict E/p require- 
ment in both decay modes to eliminate the systematic 
uncertainty from this requirement. 

The measurement of Fqoo /Fibres is based on a trigger 
that requires 25 GeV energy deposit in the Csl calorime- 
ter. As discussed in Section IIVI the minimum Ek re- 
quirement is 40 GeV, which is well above the 25 GeV 
trigger threshold. The Ek requirement ensures that 
the Kl -^ tt^tt'^tt^ acceptance is insensitive to the trigger 
threshold, but is not sufficient for K^s decays because 
of the missing neutrino. For the Ke3 acceptance to be 
insensitive to the trigger threshold, we require that the 
electron momentum be above 34 GeV/c. 

Although the Kl -^ 7r+7r~7r° decay mode is selected 
without reconstructing the tt'^ in the Csl calorimeter, the 
tt'^ decay products can hit the veto detectors. To elim- 
inate the uncertainty in modeling the veto system effi- 
ciency, F^ o/FkeS is measured using a trigger that does 

not include the veto system. Photon clusters from the 7r° 
decay may also overlap a pion cluster resulting in an E/p 
measurement that is too high. To reduce the influence 
of this effect, pion candidates are allowed to have either 
E/p less than 0.92 or E/p greater than 1.05. Allowing pi- 
ons with E/p > 1.05 recovers 1.5% of the Kl -^ tt+tt^tt^ 
sample. 
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FIG. 4; m^^ distributions after all analysis requirements ex- 
cept Tnvrvr. Selection requirements are for (a) Kl -^ tv e^u, 
(b) Kl -^ 7r±^T^ (g) Xl -^ Tr+Tr-Tr", (d) ifi -^ tt+tt". Data 
(MC) are shown by dots (histogram). "MC total" refers 
to signal plus background. "MC bkg" is the scatter- 
ing+background prediction based on simulating the other 
three (non-signal) charged decay modes. The horizontal ar- 
rows indicate the region(s) selected by the m^^ requirement. 
Note that the m,rir scale is different in (d). 



FIG. 5: pt distributions after all analysis requirements ex- 
cept pt. Selection requirements are for (a) Kl -^ tt e^u, (b) 
Kl -^ -K'^fi^v (c) Kl -^ 7r+7r"7r°, (d) Kl -^ tt+tt". Data 
(MC) are shown by dots (histogram). "MC total" refers 
to signal plus background. "MC bkg" is the scatter- 
ing-|-background prediction based on simulating the other 
three (non-signal) charged decay modes. The horizontal ar- 
rows indicate the region selected by the p^ requirement. Note 
that the pt scale is different in (d). 



B. Charged Decay Mode Background 

Table HTl summarizes the background sources that are 
subtracted from each charged decay mode. The back- 
ground to Kes decays depends on the specific selection for 
the partial width ratio, but is always less than 3 x 10~^. 
The background in the other three charged decays is 
~ 10^'^. For the partially reconstructed decay modes 
{Kl -^ TT^e^i^, Kl — » ii^n^v, and Kl -^ 7r+7r^7r°), the 

background level is checked using the k^ g distributions 

(FigureE)); based on the agreement of the data and Monte 
Carlo distributions, we assign a 20% systematic uncer- 
tainty to the background levels. The Kl -^ tt'^tt' back- 
ground is evaluated in the same manner as in the e'/e 
analysis |45| . 

There are also events in which the parent kaon has scat- 
tered in the defining collimator. This "collimator scatter- 
ing" contribution includes a regenerated iiTs-coniponent, 
and therefore must be subtracted in the Kl —^ 7r"'"7r~ 
analysis; collimator scattering is suppressed to 0.01% us- 
ing the p^ requirement. For the partially reconstructed 
decay modes, collimator scattering is suppressed to 0.1% 
using the RING requirement; this scattering component is 
not subtracted, and is included in the samples for both 
data and MC. 



TABLE II: Charged decay backgrounds from other kaon de- 
cays. 







Background 


(xlO'') to: 




Decay Mode 


K,3 


Kij,3 


K+-0 


K+- 


K,3 


— 


3 


2.4 


10 


K^3 


0.02 


— 


2.7 


5 


K^ — 


<0.1 


5 


— 


— 


K+- 


<0.2 


3 


— 


^a 


Total 


<0.3 


11 


5 


16 



"This background is from collimator scattering. 



VI. NEUTRAL DECAY MODE ANALYSIS 



Neutral Decay Mode Reconstruction and Event 
Selection 



The reconstruction of the Kl 



Ott" and 

Kl — > tt^tt^tt" decay modes is based on energies and 
positions of photons measured in the Csl electromag- 
netic calorimeter as described in [l|. Exactly four (six) 
clusters, each with a transverse profile consistent with 
a photon, are required for Kl -^ 7r°7r° {Kl -^ 7r°7r°7r°). 
The clusters must be separated from each other by at 
least 7.5 cm and have energy greater than 3 GeV. The 



fiducial volume is defined by cluster positions measured 
in the calorimeter. We reject events in which any cluster 
position is reconstructed in the layer of crystals adjacent 
to the beam holes (Fig.^fa)) or in the outermost layer 
of crystals. 

The center-of-energy of photon clusters is required to 
lie within an 11x11 cm^ square centered on the beam 
profile (RING < 121 cm^); the RING distribution for each 
neutral decay mode is shown in Fig. [7| The RING cut 
removes most events in which the kaon has scattered in 
the collimator or regenerator. 

Photons are paired to reconstruct two or three neutral 
pions consistent with a single decay vertex. The number 
of possible photon pairings is 3 for Kl -^ tt'^tt'^ and 15 
for Kl — > 7r°7r°7r°. To select the best pairing, we intro- 
duce a "pairing-x^" variable (x^o), which quantifies the 
consistency of the 7r° vertices. The pairing x^ is required 
to be less than 50 for Kl -^ 7r°7r° and less than 75 for 
Kl -^ TT^TT^TT^ |4g. This procedure identifies the correct 
photon pairing in more than 99% of the events. The Zk 
location of the kaon decay vertex is determined from a 
weighted average of the tt'^ vertices. The main kinematic 
requirement is that the invariant mass of the 2^^^ or Stt*^ 
final state (Fig.EJ) be between 0.490 and 0.505 GeV/c^. 
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FIG. 6: fc-i distributions after all analysis requirements ex- 
cept k^ 0- Selection requirements are for (a) Kl — > n e^v, 

(b) Kl -^ Tx^^'^v, (c) Kl -^ 7r+7r"7r°. Data (MC) are shown 
by dots (histogram). "MC total" refers to signal plus back- 
ground. "MC bkg" is the scattering-|-background prediction 
based on simulating the other three (non-signal) charged de- 
cay modes. The horizontal arrows indicate the region selected 
by the fc-| o requirement. 



B. Neutral Decay Mode Backgrounds 

The background subtraction procedure for Kl -^ tt^tt'^ 
is identical to that used in the e'/e analysis. The back- 
ground composition is 0.30% from scattering in the re- 
generator, 0.09% from colhmator scattering, and 0.32% 
from Kl -^ tf^tf^tt^ in which two photons are unde- 
tected ( "37r°-background" ) ji^. The total background 
is (0.71 ±0.06)%. 

In Fig. 13a), events with RING > 200 cm^ are almost 
entirely due to scattering in the regenerator and collima- 
tor; the MC predicts both the absolute level and RING- 
shape. Note that events with scattered kaons have the 
same invariant-mass distribution as events with unscat- 
tered kaons, and therefore cannot be identified in the 
7r°7r° mass distribution (Fig.lHl^a)). 

In Fig. ISIa), 97% of the events outside the signal 
region result from misreconstructed Kl — > Stt*^ events 
and the remaining 3% are from from Kl — > tt^tt'' events 
with the wrong photon pairing. These mass sidebands 
are well modeled in the simulation. The 37r*'-background 
is responsible for the apparent increase in background 
under the RING-signal in Fig.lJfa). 

In the Kl -^ Tr^Tr^Tr" decay mode, no source of back- 
ground has been identified. The contribution from kaon 
scattering (0.1% after the RING cut) is not subtracted; it 
is well modeled in the simulation, as shown in Fig. [Tfb) . 
The mass side-bands in Fig. |HIb) and (c) result from 
Kl — *• Tr'^Tr'^Tr'^ decays in which the wrong photon pair- 
ing is used to compute the invariant mass. These mis- 
pairings, which are well modeled by the MC, are not 
subtracted. 



VII. MONTE CARLO SIMULATION AND 
SYSTEMATIC UNCERTAINTIES 

A detailed Monte Carlo simulation is used to determine 
the acceptance for each decay mode. These acceptances 
are used to correct the background-subtracted numbers 
of events for each partial width ratio. The acceptance, 
Ai for decay mode "i" , is defined as 



A, = Nr^Nf 



(8) 



where Nf^^ is the number of events generated within the 
nominal Ek and Zk ranges, and TVf'^'^ is the number of 
reconstructed events 48] . To account for radiative effects 
and resolution, which can cause reconstructed events to 
migrate across the Ek and Zk boundaries, N^'^'^ is deter- 
mined from a Monte Carlo generated with broader Ek 
and Zk ranges. 

The systematic uncertainties for this analysis, which 
are summarized in Table UTTl fall into the following cate- 
gories: acceptance, background, external branching ra- 
tios, and MC statistics. Uncertainties in acceptance, 
which result from imperfections in the Monte Carlo sim- 
ulation, are by far the most important. 



TABLE III 


Systematic 


uncertainties in partial 


width ratios (in pe 


rcent) 




Source of uncertainty 


Ta'^:?/ PA-eS 


Pooo/ Pa'gS 


P+-o/ Pa'cS 


P+-/ FAe:? 


Foo/ Fooo 


Acceptance (MC Simulation) 












Event Generation: 












- Kaon energy spectrum 


0.02 


0.16 


0.04 


0.02 


0.01 


- Form factor 


0.11 


0.08 


0.29 


0.08 


0.00 


Radiative corrections: 












- 


0.15 


0.20 


0.14 


0.14 


0.00 


Particle Propagation: 












- Detector material 


0.10 


0.56 


0.33 


0.33 


0.15 


- Detector geometry 


0.02 


0.39 


0.05 


0.02 


0.08 


Detector Response: 












- Accidental activity 


0.00 


0.22 


0.04 


0.02 


0.03 


- Trigger 


0.00 


0.07 


0.10 


0.07 


0.28 


- e*,/i*,7r^ reconstruction 


0.21 


0.70 


0.24 


0.26 


0.00 


- n" reconstruction 


0.00 


0.37 


0.00 


0.00 


0.23 


Background 


0.10 


0.00 


0.02 


0.04 


0.04 


BiTv" ^ 77) 


0.00 


0.10 


0.10 


0.00 


0.03 


Monte Carlo Statistics 


0.10 


0.12 


0.05 


0.13 


0.16 


Total 


0.33 


1.12 


0.55 


0.47 


0.44 



The Monte Carlo simulation includes four main steps, 
each of which introduces systematic uncertainties in 
the acceptance: event generation, radiative corrections, 
propagation of particles through the detector, and de- 
tailed simulation of detector response. In the following 



sections, we will discuss these MC steps as well as the 
associated systematic uncertainties. 
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FIG. 7: RING distribution for (a) Kl 
Kl -^ tt^tt^tt'' candidates. Data are shown with dots. "MC 
total" (histogram) refers to the simulation of the signal includ- 
ing backgrounds and scattering in the collimator and regen- 
erator. The dashed histogram shows scattering+background 
predicted by MC. The arrow indicates the analysis require- 
ment. 



(b) FIG. 8: (a) 2^" mass, (b) 37r° mass, and (c) 3^" mass shown 



with extended mass scale. Data are shown with dots. "MC 
total" (histogram) refers to the simulation of the signal in- 
cluding backgrounds scattering in the collimator and regen- 
erator. The dashed histogram shows scattering-f background 
predicted by MC. The arrows indicate the analysis require- 
ment. 
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FIG. 9: Comparison of the vacuum beam kaon energy distri- 
butions for data (dots) and MC (histogram). For the semilep- 

tonic and K^ rnodes, the highest Ek solution is plotted. 

The data-to-MC ratios on the right are fit to a line, and the 
i?x-slopes are shown. 



A. Event Generation 

Event generation starts by selecting a kaon energy 
from a spectrum tuned with 10 million K^ — > tt+tt" 
events [l| . Figure El shows data-to-MC comparisons of 
the reconstructed energy for the four main Kl decay 
modes. To limit possible acceptance biases in the deter- 
mination of this spectrum, we compare the high-statistics 
Kl -^ vr+vr^and Kl — ^ tt^tt" spectra from the e'/e sam- 
ples. These spectra agree to better than 1% for kaon 
energies between 40 and 120 GeV (or 0.01% per GeV). 
For each partial width ratio, the systematic uncertainty 
is based on this 0.01%/GeV EK-slope uncertainty and 
the difference between the average kaon energy of the 
two decay modes. 

The decay kinematics for the 3-body decay modes 
depends on form factors. For the Kl -^ 7r+7r~7r'' and 
Kl -^ TT^TT^Tr" decay modes, PDG [3| values of the form 
factors with 2 sigma uncertainties are used. For semilep- 
tonic form factors we use our own measurements Ij- 
The largest uncertainty from form factors is in the 
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FIG. 10: Comparison of the energy distribution for radiative 
photon candidates for each of the charged decay modes. The 
data are shown as dots and the MC as a histogram. 



B. Radiative Corrections 

For Ke3 and K^^, inner bremsstrahlung (IB) contri- 
butions are accounted for using a new program, KLDR, 
described in Q- PHOTOS is used to generate the IB 
contribution in Kl —^ 7r"'"7r~7r'^ decays. The simulation 
of the Kl -^ tt+tt^ decay mode includes IB contribu- 
tions, but does not include the direct emission compo- 
nent, which has a negligible impact on this analysis. 

To check the simulation of IB, we have performed 
an analysis for each charged decay mode in which 
high-energy radiated photons are identified in the Gsl 
calorimeter. Figure ^1 shows the data and MC energy 
distribution for radiative photon candidates in the labo- 
ratory frame; the shape and normalization agree for all 
decay modes. 

The simulation without IB changes the acceptance by 
2-3% for the K^s mode, depending on the electron energy 
requirement, and less than 0.5% for the other modes. The 
systematic uncertainty in the partial width ratios is taken 
to be 6% of the acceptance change from IB based on our 
study of radiative decays. The resulting uncertainties 
vary from 0.14% to 0.20%. 



C. Particle propagation 

Once a kaon decay is generated, the decay products 
and their secondaries are propagated through the detec- 
tor. This propagation includes both electromagnetic and 
hadronic interactions of particles with detector material. 
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and requires precise modeling of the detector geometry 
and composition. To model interactions in the detector, 
GEANT is used to generate process-specific libraries 
that are used by our MC. 



1. Detector Material 

a. Electromagnetic interactions in detector material 
include photon conversions, Bremsstrahlung, multiple 
scattering, and S-ray production. For photon conver- 
sions and Bremsstrahlung, most particle losses in the 
reconstruction result from interactions upstream of the 
analysis magnet. For multiple scattering and S-iay pro- 
duction, interactions up to the last drift chamber (DC4) 
are important. We estimate 0.73% radiation lengths of 
material upstream of the analysis magnet, and a total of 
1.18% radiation lengths of material through DC4. 

To check our estimate of detector material upstream 
of the analysis magnet, we study K^ -^ 'K'^e^v decays 
in which an external Bremsstrahlung photon is iden- 
tified in the Csl calorimeter. In this study, we take 
advantage of the magnet kick to separate the electron 
from the Bremsstrahlung photon. Figure lllf a'! illus- 
trates a ii'e3 decay in which a photon is produced in 
DC2. Figure [TTT b') shows the distribution of distances 
(Ai?^brcm) between the candidate photon cluster and 
the extrapolation of the electron trajectory (measured 
in DCH-DC2) to the Csl calorimeter. The peak in 
Fig. Illf b) is mainly from external Bremsstrahlung, while 
the high-side shoulder is from radiative Ki^ — > tt e^!/7 
decays. To isolate events with external Bremsstrahlung, 
we require Ai?-),brem < 1 cm; the background from ra- 
diative Kl — > 7r^e^;^7 events is 43%. The MC sample 
is normalized to the total number of K^j^ decays in data; 
the fraction of events with a Bremsstrahlung photon in 
this study is the same in data and MC to within 5%. 

The material downstream of the analysis magnet (in- 
cluding the 0.027 Xq trigger hodoscope) is checked with 



Kl -^ Tr^Tr^TT*^ decays in which one or more of the pho- 
tons converts and gives hits in the hodoscope [i^. The 
fraction of events with hodoscope hits is measured to be 
(13.06 ± 0.03)% in data and (12.91 ± 0.03)% in MC. 

The simulation of (5-rays and multiple scattering each 
use a GEANT-based library. The J-ray simulation is tuned 
to data using DC signals recorded at earlier times than 
expected based on the track location in the DC cell. To 
check the 5-ray simulation, we consider the case in which 
a (5-ray drifts from one drift chamber cell into a neigh- 
boring cell and leaves an extra hit adjacent to the track. 
This effect is observed in 35% of the events in data (cor- 
responding to about 2% probability per wire-plane), and 
30% of the MC events. We therefore assign a 15% un- 
certainty to the (5-ray simulation. The uncertainty in our 
simulation of multiple scattering is estimated to be 10% 
of its effect on the acceptance, based on data-MC com- 
parisons of the matching of charged particle trajectories 
at the decay vertex. 
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FIG. 11: (a) Illustration of electron Bremsstrahlung at DC2 
for a Kl -^ -K e^v decay. The dotted line is extrapolated 
from the upstream electron trajectory to the Csl. The track 
bend between DC2 and DCS is from the magnet kick, (b) 
Measured distribution of A_R.^brom- Data are shown by dots, 
and MC Kei by histogram. The MC contribution from radia- 
tive Kl — > Ti^e^vy is shown by the dashed histogram. The 
MC sample is normalized to the total number of K^z candi- 
dates in data. 



b. Hadronic interactions of charged pions can re- 
sult in broad showers that lead to tracking losses and 
energy deposits in veto detectors. Pion interactions in 
the spectrometer and trigger hodoscope are simulated 
with GEANT-based libraries. Interactions in material up 
through the last drift chamber (0.7% pion interaction 
lengths) usually result in track loss. We assign a 50% 
uncertainty on this source of track loss, as described in 
Appendix ^ The corresponding 0.35% uncertainty in 
pion track-loss affects partial width ratios with a differ- 
ent number of charged pions in the numerator and de- 
nominator. 

Pion interactions in the trigger hodoscope (1.2% pion 
interaction lengths) downstream of the charged spec- 
trometer often prevent the track from matching a clus- 
ter in the Csl calorimeter. These interactions also can 
produce hadronic showers that deposit energy in the Csl- 
Veto. The inefficiency of the track-cluster match require- 
ment for pions is 0.6% for data and 0.5% for MC. The 
associated systematic uncertainties in the partial width 
ratios are negligible because only one of the two tracks is 
required to match a cluster. The Csl- Veto is only used in 
the trigger for F+./Fj^-es (Table QJ; the loss is measured 
to be (0.4 ± 0.2)% and is included in the MC. 

Most charged pions that do not decay interact hadron- 
ically in the Csl or muon system steel. The fraction of 
pions that penetrate the steel and produce a signal in 
the muon hodoscope is measured with fully reconstructed 
Kl — > 7r+7r~7r'' decays from the low intensity sample; the 
fraction is determined to be (1.0 ± 0.1) x I0~'^p.^, where 
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FIG. 12: (a) layout of inner Csl region near beam-holes. Pho- 
tons with reconstructed Csl position in a crystal adjacent to 
a beam hole (marked with an "X") are rejected in the analy- 
sis. Xb is the coordinate at one of the crystal boundaries that 
defines the photon acceptance, (b) For Kl — > tt e^u decays, 
X — Xb for electrons at the Csl determined by extrapolating 
the trajectory measured in the spectrometer. The Csl posi- 
tion requirement in (a) is applied to this electron sample. The 
200 /xm wide arrow indicates the systematic uncertainty for 
this photon aperture. 
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FIG. 13: Comparison of the vacuum beam z distributions for 
data (dots) and MC (histogram). The data-to-MC ratios on 
the right are fit to a line, and the ^-slopes are shown. 



2. Detector Geometry 

The dimensions of the four drift chambers are known 
to better than 20/ini based on optical surveys. The spec- 
trometer is aligned in situ as explained in 1] . The Csl in- 
ner aperture for photons is illustrated in Fig. ll2r a). This 
aperture is measured by comparing extrapolated electron 
tracks measured in the spectrometer with cluster posi- 
tions measured in the calorimeter (Fig. ITW b)): the un- 
certainty in this aperture size is 200/im. The calorimeter 
dimensions (1.9 m x 1.9 m) are known to better than 
1mm. 

The agreement of data and MC decay vertex distri- 
butions (Fig. I13|l provides a sensitive overall check of 
the detector geometry. For each partial width ratio, the 
systematic uncertainty is the product of the data/MC 
Zx-slope and the difference between the average recon- 
structed Zk of the two decay modes. 



D. Detector Response 

The Monte Carlo includes a detailed simulation of de- 
tector response to different particle species, as well as the 
effect of accidental activity. 



1. Accidentals 

Accidental detector activity is measured with a trigger 
that is proportional to the instantaneous beam intensity. 
We overlay an accidental data event on each generated 
decay in the Monte Carlo simulation. Based on data- 
MC comparisons of extra detector activity, the system- 
atic uncertainty is estimated to be 10% of the acceptance 
change arising from overlaying accidental data events in 
the simulation. The effects of accidentals result in 0.22% 
uncertainty in the TaQo/TKeS ratio, and less than 0.04% 
in the other partial width ratios. 



2. Trigger 

As described in [J, the KTeV MC includes a detailed 
trigger simulation. Although the analysis requirements 
are designed to be stricter than the trigger, potential 
anomalies in the trigger pre-scales and signal process- 
ing could result in losses that are not simulated. The 
general strategy to estimate these effects is to use K^s 
decays from a charged-track trigger to study the total 
Csl energy trigger (Table P) , and to use charged decay 
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modes collected in the total-energy trigger to study the 
charged-track trigger. 

For partial width ratios in which events for the nu- 
merator and denominator are collected with the same 
trigger, uncertainties in trigger efhciency largely can- 
cel, and resulting uncertainties are less than 0.1%. For 
Foo/Foooj the only ratio measured with two separate trig- 
gers, the systematic uncertainty from the trigger effi- 
ciency is 0.28%. 



3. Response to Charged Particles 

Each charged decay mode includes a different combi- 
nation of e , TT , and /i . The simulation of the drift 
chamber response includes measured wire inefficiencies, 
and several subtle effects that cause non-Gaussian tails 
in the position resolution [1| . To simulate the response in 
the Csl calorimeter, a separate GEANT library is gener- 
ated for e^ , vr^ , and fi^ . Tails in the Csl energy response 
are measured in data as described in Appendix IbI 

Sources of systematic uncertainty are the tracking ef- 
ficiency, drift chamber calibration, tails in the E/p dis- 
tribution, and analysis cut variations. The uncertainty 
in the drift chamber efficiency is 0.6% as explained in 
Appendix^ this uncertainty affects only the Tqoq /TKeS 
ratio. The DC calibration introduces a systematic uncer- 
tainty less than 0.1% on each ratio. The effect of tails in 
the Csl energy response introduces systematic uncertain- 
ties well below 0.1% on the charged partial width ratios. 
Cut variation studies introduce a 0.2% uncertainty on all 
partial width ratios. 



4- Response to Photons 

The Csl calorimeter response to photons is simulated 
with a GEANT library. The low-side tail in the energy 
response is assumed to be the same as for electrons (Ap- 
pendix . The most crucial role of the photon simu- 
lation is to predict the efficiency of reconstructing the 
Kl -^ n^n^TT^ decay mode for the Tooo/^Ke3 ratio. The 
sources of systematic uncertainty in reconstruction of 
multi-TT*' events are: photon pairing efficiency, energy 
scale, and photon reconstruction efficiency. 

The pairing efficiency study uses Kl -^ Tr^Tr^Tr" decays 
since there is no background after identifying six pho- 
ton clusters in the Csl calorimeter. The sidebands in 
the Stt" invariant mass distribution, as well as events 
with x^a values beyond the selection cut, result from 
misreconstructed Kl —^ tt^tt'^tt'^ events. These misrecon- 
structed events result almost entirely from selecting the 
incorrect photon pairing. Fig. |HIc) illustrates misrecon- 
structed events in the Stt*' invariant mass distribution. 
The MC sample in Fig. |SJc) is normalized to data in the 
15 MeV wide signal region; the data and MC sidebands 
are in good agreement. The combined requirements on 
X^o and 37r° invariant mass remove (0.80 ±0.01)% of the 



Kl -^ TT^TT^TT^ events in data, and remove (0.66 ±0.01)% 
of the events in MC. This 0.14% difference is included as 
a systematic uncertainty on the Fooo/FA'eS ratio. 

Uncertainties in the calorimeter energy scale and lin- 
earity affect the 7r° reconstruction efficiency, primarily 
because of the photon energy and kaon energy require- 
ments. The systematic uncertainty is based on the same 
set of tests as in the e'/e analysis; these tests result in a 
0.33% uncertainty on the Fooo/Fifes ratio, and a 0.05% 
uncertainty on the Foo/Fqoo ratio. 

We consider three sources of uncertainty in the photon 
reconstruction efficiency: detector readout, dead mate- 
rial, and Csl cluster shape requirements. The calorime- 
ter readout inefficiency is monitored with a laser system, 
and is measured to be less than 10~^. The amount of 
dead material between crystals is checked with muons, 
and is included in the MC; the probability of losing a 
photon because of this dead material is less than 10~^. 
The effect of the photon cluster shape requirement for 7 
candidates in the calorimeter is studied by removing this 
cut in the Fooo/Fa-eS analysis; the corresponding change 
of 0.05% is taken as a systematic uncertainty. 



E. Sensitivity to vr" Branching Fractions 

For decay modes that use tt" -^ 77 decays, we correct 
for the branching ratio, B(7r° -^ 77) = 0.9880 ± 0.0003. 
The decay 7r° — > e"'"e~7 has a negligible effect on all decay 
mode acceptances except for K^ -^ 7r+7r~7r°. Although 
the measurement of Kl -^ tt^tt^tt'^ ignores the tt", it is 
still sensitive to the tt" — > e+e~7 decay because the extra 
tracks can cause the event to be rejected. This effect is 
studied using a MC sample of Kl -^ n^n^n^ decays in 
which the 7r° decays to e~^e~j; the acceptance for this 

MC sample is 0.32 x A^ Oi where A^ is the nominal 

acceptance for Kl -^ tt^tt^tt^ with 7r° -^ 77. We assign 
a 20% uncertainty to the fraction of these events pass- 
ing the selection, resulting in a 0.1% uncertainty in the 
F^ o/FxeS partial width ratio. 



VIII. RESULTS 
A. Partial width ratios 

The numbers of events and detector acceptances for 
all decay modes are summarized in Table HVI the result- 
ing partial width ratios are given in the last column of 
the table. For each partial width ratio, the first error 
is statistical and the second systematic. The systematic 
uncertainty is calculated as the sum in quadrature of the 
individual sources. Note that although the partial width 
ratios use independent data samples, there are correla- 
tions among the systematic errors. For example, uncer- 
tainties from external bremsstrahlung cause correlated 
uncertainties among the four partial width ratios involv- 
ing TkcS- These correlations are treated as described in 
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Appendix D of Ref. [l|. The correlation coefficients for 
the five partial width ratios are given in Table IVI 

The systematic precision of each ratio depends on the 
cancellation of systematic uncertainties for the pair of 
modes. Among the partial width ratios, the measurement 
of rooo/r_R-e3 has the largest systematic uncertainty, with 
Csyst — 1.12%. This ratio has the largest uncertainty be- 
cause there is no cancellation of the 0.6% uncertainty in 
the tracking efhciency, as well as the 0.37% uncertainty 
in the Stt" reconstruction. Pairs of modes with the same 
number of charged pions in the final state {^Kfis/^KeS 
and Foo/rooo) have the smallest systematic uncertainty, 
with a 



syst 



0.4%. Pairs of modes in which the number 

of charged pions is different (T^ q/^KcS and P^i /^Kes) 

have a larger systematic uncertainty of ffsyst '^ 0.55%; 
this increased uncertainty is mainly from the 0.35% un- 
certainty in losses from hadronic interactions. 



B. Cross Checks of Partial Width Ratios 

We have performed several crosschecks of our measure- 
ments. Some of these checks affect all of the partial width 
ratios, while others are relevant to a specific decay mode. 

As discussed in Sec. lIIII we measure each partial width 
ratio in both high and low intensity data samples. The 
factor of 10 difference in beam intensity between these 
samples results in significantly different tracking and pho- 
ton cluster reconstruction efficiencies. For example, the 
tracking inefficiency is nine times smaller in the low in- 
tensity sample (0.38% vs. 3.3%). FigureirHcompares the 
partial width ratios measured in these two samples. The 
measurements are in good agreement: the x^ P^r degree 
of freedom is 3.4/5. 



• high intensity beam with regenerator 
O xlO lower intensity without regenerator 



('"oo'^ooo' h 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 

decay rate ratio (arbitrary scale) 

FIG. 14: Partial width ratios measured with high intensity 
(solid dots) and with low intensity (open circles). Each ratio 
is normalized so that the quoted result is one. The error bars 
reflect the statistical uncertainties between the two samples. 

To check our nominal Kj^ — > tt /i^z/ analysis, which 
does not use the muon system, we perform an anal- 



ysis that requires one track to be matched to hits in 
the most downstream muon hodoscope. Since the muon 
track is identified, there are only 2 kaon energy solu- 
tions instead of the 4 solutions in the nominal analysis. 
This alternate analysis differs from the nominal analy- 
sis by (0.08 ± 0.02gig^i)%. It is worth mentioning that 
requiring a signal in the muon hodoscope does not sig- 
nificantly reduce the background because K^ — > Tr+Tr^Tr*^ 
and Kl —> tt+tt" decays can pass the K^^ selection only 
if one of the pions decays in flight. 



The nominal Kj^ 



"TT tt" analysis, which does not 



reconstruct the tt , is checked by performing an anal- 
ysis in which the Kl -^ tt+tt^tt'' decay is fully recon- 
structed using the tt'' -^ 77 decay. Requiring a recon- 
structed TT*^ —I- 77 decay in the Csl calorimeter reduces 
the acceptance by a factor of four. To increase the sta- 
tistical significance of this crosscheck, we use an indepen- 
dent sample (with x5 smaller pre-scale) collected in the 

trigger used to measure F_| /F^es (see TablePi; the two 

methods agree to (0.03 ± 0.28stat)%- 

The stability of the results is also tested by dividing the 
data into a variety of subsamples based on criteria such 
as vertex position, kaon energy, minimum track separa- 
tion, and minimum photon energy. The measured partial 
width ratios are found to be consistent within the uncor- 
related statistical uncertainty in all of these studies. 

The clean separation of semileptonic and 

Kl -^ Tr+vr^vr*' decays in the fc^ distribution 

(Fig. EI^)) allows a measurement of R^ — 

F^ q/(Fx/j3 + Fxe3 -t- F^ 0) that does not use any par- 
ticle identification information. The difference between 
this fit and the nominal analysis is (0.35 ± 0.51)%. 

Assuming lepton universality, we can make an inde- 
pendent prediction of the Tjifj.3/^Ke3 ratio: 



Kfj.3 



Tk 



e3 



1 



5'' 



pred 



1 



^K 



jM 
^K 



'if 



(9) 



Here, 5j^ represents the mode-dependent long-distance 
radiative correction to the total decay width, and /J^ and 
/^ are the decay phase space integrals, which depend 
on the form factors. Using the KTeV measurement of 
/^//|, = 0.6622 ± 0.0018 Bland (1 + (5^)/(l + S%^) = 
1.0058 ±0.0010 from KLORji, we find that the ratio of 
the directly measured to predicted values of Tn^is/^KeS 
is 0.9969 ± 0.0048, consistent with 1. 



C. Determination of Branching Fractions, Partial 
Decay Widths, and |r;+ | 

Imposing the constraint that the sum of the six largest 
branching fractions is 0.9993, we determine the branch- 
ing fractions shown in Table IVII Correlations among 
the partial decay width measurements (Table are 
taken into account in calculating uncertainties in the 
branching fractions. The correlation coefficients for the 
six branching fractions are given in Table IVIII Us- 
ing the PDG average for the neutral kaon lifetime |50l |. 
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TABLE IV: Background-subtracted numbers of events, detector acceptances, and resulting ratios of partial decay widths. The 
XX next to a number of events reflects a prescale that must be applied to calculate the partial width ratio. For the partial 
width ratios, the flrst error is statistical and the second systematic. Note that differences in event selection requirements for 
the partial width ratios result in a range of acceptances for Kl -^ tt^e^u. 



Decay Modes 



Numbers of Events 



Acceptance 



Partial Width Ratio 



rooo/rife3 

r+-o/rKe3 
r+-/rife3 

roo/Eooo 



394300/449379 


0.239/ 0.180 


209365/211950 


0.046/ 0.022 


799501/(807343 x 2) 


0.200/ 0.124 


83725/(979799 x 8) 


0.265/ 0.121 


100365/(1609324 x 5) 


0.150/ 0.054 



0.6640 ± 0.0014 ± 0.0022 

0.4782 ± 0.0014 ± 0.0053 

0.3078 ± 0.0005 ± 0.0017 

(4.856 ± 0.017 ± 0.023) x 10"=* 

(4.446 ± 0.016 ± 0.019) x 10"^ 



TABLE V: Total uncertainties and correlation coefficients for the partial width ratios. 




rA>3/ ^Ke3 Looo/ LxeS ^+-o/ T^KeS ^ + -/ ^KeS 


Loo/ Looo 


Total Error 


0.0026 0.0055 0.0018 0.029 x lO"'' 


0.025 x 10^^ 



r/f^3/ rKe3 


1.00 


Looo/ LifeS 


0.14 


Fh — o/ FifeS 


0.21 


F + -/ FKe3 


0.24 


Foo/ Fooo 


0.09 



1.00 

-0.06 
-0.07 
0.30 



Correlation coefficients 



1.00 
0.49 
0.04 



1.00 
0.07 



1.00 



TABLE VI: Kl branching fractions and partial widths (Fi). 
The partial width measurements use the PDG average for the 
Kl lifetime: tl = (5.15 ± 0.04) x 10"® sec 0. The quoted 
errors are the sum in quadrature of statistical and systematic 
uncertainties. 



Decay Mode Branching Fraction 



F, (lO^g-^) 



Kl- 


-+ TT^e^U 




0.4067 ±0.0011 


Kl- 


-* n jj.^i' 




0.2701 ± 0.0009 


Kl- 


-. TT+^-V' 




0.1252 ±0.0007 


Kl- 


-. ttV^o 




0.1945 ±0.0018 


Kl- 


-^ tt'^^tt" 


(1.975 ±0.012) x 10-3 


Kl- 


-.tt^ttO 


(0, 


.865 ±0.010) x 10-3 



0.7897 ± 0.0065 
0.5244 ± 0.0044 
0.2431 ± 0.0023 
0.3777 ± 0.0045 
(3.835 ±0.038) x IQ-^ 
(1.679 ±0.024) X 10^3 



Tl — (5.15 ± 0.04) X 10 ® sec, these branching fractions 
correspond to the partial widths quoted in the same ta- 
ble. 

The Kl — > TTTT measurements, combined with the kaon 

lifetimes, also provide a precise measurement of jy^-i p = 

r{KL -^ TT+7r-)/T{Ks -^ TT+TT-y. 



TS 



B' 



B' 



[I ± &Re{e'/e)] 



TL 



B^iv 



(10) 



where B^^ _ and B\ o are the Kl — ^ tttt branching 



fractions quoted in Table |Vl| rs = (0.8963 ± 0.0005)x 
10-1° sec jHH, and i?e(e'/e) = (16.7±2.3) x 10-^ QllQ- 
We calculate the Ks -^ nlv branching fraction B^g^ = 
0.118% assuming that r{Ks -^ niv) = T{Kl ~^ ttIi^). 
The resulting value of \ri^ | is 



w^ 



= (2.228 ±0.010) X 10" 



(11) 



The uncertainty in tl contributes 0.009 x 10 '^ to the un- 
certainty in 1?;^ I while our Kl — > tttt branching fraction 



measurements contribute an uncertainty of 0.005 x 10 ^; 
Ts and i?e(e'/e) contribute negligibly to the error in jry-j |. 



IX. COMPARISON WITH PREVIOUS PARTIAL 
WIDTH AND |?7+_i MEASUREMENTS 

The new KTeV measurements of the partial width ra- 
tios and Kl branching fractions are on average a factor 
of two more precise than the current world average val- 
ues, but are not in good agreement with these averages. 
FigureElshows a comparison of the KTeV and PDG val- 
ues for the five partial width ratios. Of the five partial 
width ratios, only the Foo /Fooo measurement is in good 
agreement; note that Foo/rooo is the only ratio that does 
not include the Kl — > tt e^iy decay mode. Our measure- 
ments of TKfis/T^KeS and r_| o/TxeS disagree with the 

PDG by 5%, and our measurements of Tooq/TkcS and 
r+_/rxe3 disagree with the PDG by 10%. Figure [H 
shows the corresponding comparison of KTeV and PDG 
branching fractions. The discrepancies between KTeV 
and the PDG can be reduced significantly by applying a 
7% relative shift to either the KTeV or PDG values for 



B{Kl 



7r±e^ 



^)- 



Another measurement, which has not yet been in- 
cluded in the PDG summary, is the KLOE measurement 
of i?5 = T{Ks -^ TT+TT-)/r{Ks -> ttOttO) = 2.236 ± 
0.015 11]. The KTeV measurements of B{Kl -^ tt+tt-) 
and B{Kl — > tt^tt^), along with the world average value 
of i?e(e'/e), give Kg = 2.261 ± 0.033 in good agreement 
with, but less precise than, the KLOE result. 

To understand the discrepancy between KTeV and the 
PDG averages, we have considered the 49 measurements 
and fit results used in the PDG averages Ql- 34 of 
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TABLE VII: Total uncertainties and correlation coefficients for the Kl branching ratios. 



Total Error 



0.0011 



0.0009 



0.0018 



0.0007 



0.012 X 10" 



0.010 X 10" 



Correlation coefficients 



B(Kl 


^ 7r±eT/.) 


1.00 










B{Kl 


^7r±/iTz/) 


0.15 


1.00 








B(Kl 


-^ 7r%%«) 


-0.77 


-0.62 


1.00 






B(Kl 


-^TT + ^-TT^) 


0.18 


0.08 


-0.54 


1.00 




B(Kl 


-^TT + TV-) 


0.28 


0.22 


-0.48 


0.49 


1.00 


B{Kl 


-^ TT TV ) 


-0.72 


-0.54 


0.89 


-0.46 


-0.39 



1.00 



these measurements involve decay modes with branch- 
ing fractions greater than 1%. Figure [T7I shows the dis- 
tribution of residuals (normalized by uncertainty) be- 
tween these measurements and values obtained using the 
new KTeV measurements. The x^ P^r degree of free- 
dom is 82.9/34 showing a clear inconsistency. Approxi- 
mately 40 units of the x^ come from the three measure- 
ments with greater than 3<7 disagreement with KTeV: 
the measurement of TKfis/^KeS reported by Cho 80 Q, 

and the measurements of Fooo and r_| o/^KeS from 

NA31 10]. It is interesting to note that in the same pa- 
per, NA31 reports rooo/r+-o = 1.611±0.037, not involv- 
ing the Kl -^ -K^e^v decay mode, which is consistent 
with KTeV's measurement of rooo/r+-o = 1.567±0.020; 
this NA31 measurement is not used in the PDG branch- 
ing ratio fit because it is not independent of the other 
NA31 measurements included in the fit. 

Figure [TSl compares our new determination of |ry^ | 

with the two measurements based on K^-Ks interference 
[53 ■ The average of these two previous measurements 

gives |?7-| I = (2.295 ± 0.025) x 10~^, which disagrees 

with the KTeV evaluation by 2.7tT. Figure^|also shows 

\rj^ I determined from the charge asymmetry assuming 

CPT invariance [53 ; the value is consistent with all other 
measurements. 
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FIG. 15: Partial width ratios measured by KTeV (dots) and 
from PDG fit (open circles). 



X. CONCLUSIONS 



(2.228 ±0.010}x 10-3, reported in this paper, and |Ks|, 
described in 1391. 



In summary, we have measured the branching fractions 
for the six K^ decay modes with branching fractions 
greater than 0.05%. The new measurements are about 
a factor of two more precise than current world aver- 
ages, but are not in good agreement with these averages. 
Compared to the PDG fit 2], the KTeV measurement 
of b\Kl -^ TT^e^i/) is higher by 5%, B{Kl -^ tt°tt°it°) 
is lower by 8%, B{Kl — *■ Tr+Tr") is lower by 5%, and 
B{Kl -^ TT^TT^) is lower by 8%. Our measurements of 



are consistent 



B{Kl -^ -K^ijL^v) and B{Kl 
with the PDG fit. 

The new K^ branching fractions will require the ad- 
justment of several rare K^ branching fractions for which 
the main K^ decay modes are used as normalization. 
The Kl branching fraction measurements also may be 
used to determine several parameters, including |?7^ | — 
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APPENDIX A: DETERMINATION OF 
TRACKING LOSSES 

The determination of the absolute track reconstruc- 
tion efficiency is most important for the measurement 
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FIG. 16: Kl branching fractions measured by KTeV (dots) 
and from PDG fit (open circles). 



of ^aoa/^KeS which compares a decay mode with two 
charged tracks to a mode without charged particles in the 
final state. The causes of tracking loss fall into four cate- 
gories: (i) missing or corrupted hits induced by accidental 
activity, S-iays, and non-Gaussian tails in the chamber 
response; (ii) many spurious hits, mainly from accidental 
activity, which confuse pattern recognition; (iii) failure of 
track segments to satisfy matching criteria at the anal- 
ysis magnet or decay vertex because of large angle scat- 
tering and IT — > /ii/ decays; (iv) hadronic interactions in 
the spectrometer. All of these effects are included in the 
Monte Carlo simulation, as described in [2]. Pion decay 
is the main source of track loss, and is well described in 
the MC. This appendix describes an inclusive study of 
all other sources of track loss. The first three sources of 
track loss are related to the drift chamber performance, 
and are described in Sec. lA II The effect of hadronic 
interactions is described in Sec. IA2I 

The measurement of the tracking loss is based on par- 
tially reconstructed Kl — > 7r+7r~7r° decays that are iden- 
tified by the tt" — > 77 decay along with one or two addi- 
tional hadronic clusters in the Csl calorimeter (i.e, clus- 
ters with transverse profile inconsistent with a photon). 
Ideally, Kl — > 7r+7r~7r'' decays could be cleanly identified 
using only the four clusters in the Csl calorimeter, leading 
to a direct measurement of the two-track inefficiency; un- 
fortunately, this sample has significant background. To 
obtain a Kl — *• Tr+Tr^Tr" sample with low enough back- 
ground for this study, we also require either a completely 
reconstructed pion track or some reconstructed track seg- 
ments. With this partial reconstruction, the missing 
track information can be predicted from kinematic con- 
straints, and then compared with the track information 
found by the track reconstruction. The data used for this 
study were collected in a trigger that requires energy in 
the Csl calorimeter and rejects events with activity in 
the muon system. 

The two-track inefficiency is measured in two steps. 
First, we measure the single track inefiiciency after one 
track is reconstructed (rji). Next, we measure the proba- 
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FIG. 17: Residual distribution between 34 individual mea- 
surements (Mexp) used in the PDG fit and the KTeV determi- 
nations of these quantities (Mktbv). o"am is the uncertainty 
on MktcV — Mexp. Measurements disagreeing with KTeV by 
more than 3a are shown as solid squares. The measurements 
are Guidoni 65 [H, Astbury 65 12], Hawkins 66 O, 
Hopkins 67 |3, Kuiyukina 68 ilSil, Evans 73 Ea, Bran- 
denburg 73 117], Alexander 73 Till. Cho 7701, NA31 
95 10], Aieksanyan 64 20], Budagov 68 "21], Anikina 

64 22], Williams 74 23], Cho 80 9], Franzini 65 24], 
Cho 70 [25], Webber 71 26], Burgun 72 |2l, Aubert 

65 111], Chan 71 ^, Behr 66 ^, Webber 70|3l|l. Meis- 
ner 71 ^, Cho 71 [13, James 72 ^, Baldo-Ceolin 
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FIG. 18: Comparison of KTeV's {rj^ ] measurement with 

previous measurements based on Kl-Ks interference |3a . l37t . 
and the semileptonic charge asymmetry p.l38l|. Note that the 

CPLEAR value of \r)^ ] has been adjusted to ts ~ 0.8963 x 

10"^" sec using their quoted dependence on the Ks lifetime. 



bility of finding exactly zero tracks (770)- The correspond- 
ing two-track loss, 772, is then given by 2rii + rjo. 
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1. Tracking Inefficiency 

The study of the single track efRciency uses events 
with a reconstructed tt'^ —> 77 along with two hadronic 
clusters in the Csl calorimeter. Requiring exactly two 
hadronic clusters suppresses Kl -^ 7r"'"7r~7r° decays with 
hadronic interactions in the spectrometer because such 
events tend to create additional hadronic clusters. Events 
with TT ~> fii' decays are also suppressed because these 
events leave only one hadron cluster. One of the two 
hadron clusters must be matched to a fully reconstructed 
track in the spectrometer, leaving two possible kine- 
matic solutions for the missing track. The fully re- 
constructed track, along with the measured position of 
the other hadron cluster (associated with the missing 
track), provides sufficient information to select the cor- 
rect kinematic solution for the missing pion trajectory. 
FigurelT^ illustrates this selection. For the high intensity 
sample, we find the single track inefficiency (rji) to be 
(1.47 ± 0.02)% in data and (1.32 ± 0.02)% in MC; in the 
low intensity sample, 771 = (0.18 ± 0.02)% in data and 
(0.13 ±0.02)% in MC. 

In the above study, the reconstructed track require- 
ment excludes the case in which correlated hit losses 
within a single drift chamber result in no reconstructed 
tracks. To account for correlated losses within a drift 
chamber, we perform a separate analysis to measure the 
fraction of events with zero reconstructed tracks (jiq). 
We start with the same calorimeter selection of two pho- 
ton clusters for 7r° — » 77 and two hadron clusters. Next, 
we require two reconstructed track segments, either both 
upstream of the analysis magnet in DCl and DC2, or 
both downstream of the analysis magnet in DC3 and 
DC4. Figure 1201 shows the Tr+Tr^Tr" invariant mass dis- 
tributions for these two samples. The 7r"'"7r~7r'^-mass res- 
olution for the DC1-DC2 (DC3-DC4) selection is about 
3 MeV/c^ (2 MeV/c^) compared to 1 MeV/c^ for the 
standard reconstruction. The tails are well described by 



the Kl 



'n" MC. For the high intensity sample. 



we find r]o = (0.35 ± 0.02)% in data and (0.06 ± 0.02)% 
in MC. For the low intensity sample, 7/0 = (0.03 ±0.01)% 
in data and zero in MC. 

In the analyses for both 771 and ryo, we ensure that the 
special selection does not introduce biases in the tracking 
efRciency, and that the efficiency corrections are applica- 
ble for the nominal selection of charged decay modes. 
One of the most important aspects in the tracking ineffi- 
ciency study is to require separation between the recon- 
structed track (or track segment) and the kinematically 
predicted track (or track segment). This requirement is 
necessary because nearby tracks have a large tracking 
inefficiency, which is irrelevant for the nominal analysis 
that uses the track separation cut. Figure [2] shows 771 as 
a function of X-separation (AX) between the two tracks 
at DCl. At small track separation, the inefficiency is 
larger than 20%. A requirement of AX > 7 cm safely 
removes tracks with small separation. 

The total two-track inefficiency (772 = 27yi+7yo) depends 



strongly on beam intensity. For the high intensity period, 
772 = (3.28 ± 0.04)% in data; it is not completely repro- 
duced by the simulation for which 772 = (2.70 ± 0.05)%. 
For the low intensity period, the inefficiency in data is 
much smaller, (0.38 ± 0.06)%; the MC inefficiency is 
(0.26 ± 0.05)%, which agrees with the data. Based on 
these measurements, we apply a —0.6% correction to the 
acceptance for the high intensity period and a —0.12% 
correction for the low intensity period. The uncertainty 
on this correction is taken as 100% of the effect: 0.6% 
and 0.12% for the high and low intensity periods, respec- 
tively. 



2. Hadronic Interactions in the Spectrometer 

From our Monte Carlo simulation, the track loss from 
hadronic interactions in the spectrometer is 0.7% per 
pion track. Since these interactions also cause the Csl 
cluster associated with the track to be lost, we cannot 
measure the track loss using the method described in 
Sec. lA II To identify hadronic interactions clearly with- 
out the second pion cluster, we tag pions with an in- 
teraction between DC3 and DC4. We use partially re- 
constructed Kl — > 7r+7r~7r° decays using 7r° -^ 77, one 
complete pion track, and hits from the second pion track 
in the first three drift chambers. Note that only one 
hadronic cluster, associated with the complete track, is 
required. 

This study uses low intensity data for which the track- 
ing inefficiency has a small effect. Events are vetoed if 
there are extra hits in any of the first three drift cham- 
bers; this further suppresses accidental activity and sim- 




track inefficiency {Z-X projection). Solid lines represent par- 
ticles reconstructed in the detector. The dashed line indicates 
the pion trajectory calculated by the kinematics of the recon- 
structed particles with the assumption of a Kl -^ tt^tv'tt'^ 
decay. 
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plifies the particle trajectory determination. Hadronic 
interactions are tagged by requiring no hit in the last 
drift chamber within 1 cm of the second track extrapola- 
tion; according to MC, 90% of these events are due to a 
hadronic interaction. In data, (0.14±0.02)%of events are 
tagged as having hadronic interactions; the correspond- 
ing fraction for MC is (0.096 ± 0.005)%. The difference 
between these two fractions leads to a 50% uncertainty 
on the track loss from hadronic interactions. 



APPENDIX B: TAILS IN THE CSI ENERGY 
RESPONSE 



To measure the partial width ratios, it is important to 
understand tails in the Csl energy response for the four 
different types of particles detected in the analysis: elec- 
trons, photons, pions, and muons. The energy response 
tails affect the absolute detection efficiency, as well as 
particle misidentification that leads to background. In 
this appendix, we describe these tails, how they are de- 
termined from data, and how they are simulated. 

KTeV's GEANT-based MC does not include the effects 
of photo- nuclear interactions (jN). If a jN reaction 
occurs during the electromagnetic shower development, 
a neutron or charged pion can escape the calorimeter, 
resulting in an energy deposit that is too low. Im- 
perfections in the treatment of dead material (wrap- 
ping) between crystals also can lead to energy tails 
that are not properly modeled. These effects are em- 
pirically modeled and included in the MC by fitting a 
function to the low-side tail in the electron E/p dis- 
tribution (Fig. |2Ia)). The key point in this procedure 
is to select a sample of K^ -^ ir^e^v decays with low 



|10 5 


(a) .DC1+DC2 


•* 


•A 




0) 


• data 

- MC ■ 






lS^qS 


r * 

J \^ 




A 1 "^te»i«« 




m\ 






10 2 




1 


1 




0.45 0.475 0.5 0.525 0.55 0.45 0.475 0.5 0.525 0.55 



iCn 71 mass (GeV/c 



%*% %° mass (GeV/c^ 



FIG. 20: For 770 tracking inefficiency study, 7r+7r~7r° in- 
variant mass distributions after all analysis requirements ex- 
cept 7r^7r~7r''-mass. The data are shown as dots and the 
MC as a histogram. Distributions are based on identifying 
track segments in (a) DC1-I-DC2 and (b) DC3-I-DC4. Events 
with invariant mass above 0.55 GeV/c^ are shown in the last 
bin. The horizontal arrows indicate the region selected by 
the 7r"'"7r~7r° mass requirement. The rightmost bin shows all 
events with mass greater than 0.55 GeV/c^. 



enough background to avoid pion contamination for elec- 
tron E/p values as low as 0.6. To achieve such elec- 
tron purity, the pion is required to satisfy E/p < 0.3 
to avoid swapping the electron and pion assignments, 
and TOtttt > 0.370 GeV/c^ to reject background from 
Kl -^ TT+Tr^TT*^ decays. Kl — > n^fi^iy decays are vetoed 
by the muon system. With the nominal electron selection 
requirement of E/p > 0.92, the inefficiency is measured 
to be 0.2%. This non-Gaussian tail also affects photons in 
the neutral decay modes; the presence of this tail changes 
the Kl -^ 7r°7r°7r° acceptance by 0.7%. 

The electron E/p tail in Fig. |2Ia) shows losses up to 
40% of the incident energy (i.e., down to E/p ~ 0.6). 
To check for anomalous energy losses of more than 40%, 
we use a sample of muons collected in dedicated runs as 
explained in Q ■ We select muons that are at least 3 mm 
away from any crystal boundary to avoid the known ef- 
fects of the small amount of dead material (wrapping) 
between crystals. Figurel^shows the energy deposit dis- 
tribution for 5.5 million muons hitting the Csl calorime- 
ter away from crystal boundaries. The fraction of events 
with an anomalously low energy deposit (below 200 MeV) 
is 5 X 10^^ in data; this effect is negligible in the branch- 
ing fraction analysis. 

For charged pions in the Csl calorimeter, the E/p dis- 
tributions for data and GEANT-based MC agree reason- 
able well (Fig. Efb)). The pion E/p cut inefficiency is 
measured in both data and MC with Kj^ -^ 7r"'"7r~7r° de- 
cays in the low intensity sample. Note that with the 
full reconstruction of Kl — > 7r+7r~7r°, the pion E/p cut 
is not needed to achieve a negligible background. The 
inefficiency difference between data and MC is measured 
as a function of the E/p requirement and the proximity 
to the beam holes. The data-MC inefficiency difference 
varies between 0.1% and 0.3%, and is used to correct the 
MC samples. 

In the Kl — > tt fi^i' analysis, at least one track is re- 
quired to deposit less than 2 GeV in the Csl calorimeter 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 2 
AX at DC1 (m) 

FIG. 21: Single track reconstruction inefficiency (r;i) as a 
function of the X-distance (AX) between the reconstructed 
track and the calculated track at DCl. The arrow indicates 
the selection requirement. 
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(i.e., less than 5 times the energy deposit of a minimum 
ionizing particle). In a separate K^j^^ analysis that identi- 
fies the muon with the muon hodoscope, the inefficiency 
of this cluster energy requirement is measured to be 0.3% 
in both data and MC. No correction is used to simulate 
the muon energy response in the Csl calorimeter. 
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FIG. 22: Energy deposit in the Csl for muons with track 
projection at the Csl at least 3 mm from any crystal boundary. 
The data are shown as dots and the MC as a histogram. 
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